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We model a microscopic heat engine as a particle hopping on a one- dimensional lattice in a periodic 
sawtooth potential, with or without load, assisted by the thermal kicks it gets from alternately placed 
hot and cold thermal baths. We find analytic expressions for current and rate of heat flow when 
the engine operates at steady state. Three regions are identified where the model acts either as a 
heat engine or as a refrigerator or as neither of the two. At quasistatic limit both efficiency of the 
engine and coefficient of performance of the refrigerator go to that for Carnot engine and Carnot 
refrigerator, respectively. We investigate efficiency of the engine at two operating conditions (at 
maximum power and at optimum value with respect to energy and time) and compare them with 
those of the endoreversible and Carnot engines. 

PACS numbers: 5.40. Jc Brownian motion-05.60.-k Transport procssses -05. 70. -a Thermodynamics 



I. INTRODUCTION 

Even though both macroscopic as well as microscopic 
heat engines work on the same thermodynamic princi- 
ples, wide-ranging studies have been done in improving 
the performance of macroscopic heat engines jlj. At 
present, the study of microscopic heat engines has re- 
ceived considerable attention 0,0] • This is because of the 
trend in miniaturization and the need to utilize energy 
resources available at microscopic scales. As such, mod- 
elling microscopic heat engines and finding how well they 
perform is a primary task to be undertaken at present. 

To get a first insight as to how such engines perform, 
it is important to take a toy model that has the basic 
ingredients. In a recent paper Qj we considered a simple 
model of a Brownian heat engine and found exact ana- 
lytic expressions for quantities like current, efficiency and 
coefficient of performance. This, in turn, enabled us to 
explore different features of the engine such as efficiency 
at maximum power, optimized efficiency as well as effi- 
ciency at quasistatic limit. The present work addresses 
the same basic issues of a tiny heat engine. However, 
here the particle moves on a discrete lattice by hopping 
as opposed to a continuous Brownian motion in a viscous 
medium. Even though the model and its corresponding 
dynamics is completely different from the one previously 
studied, the results can be compared with those of the 
previous work at least qualitatively. As such, this work 
can be taken as an independent check of the results found 
in the previous work on microscopic heat engine. 

The paper is organized as follows: In section II, we will 
first introduce our model in the absence of external load 
and set up the dynamics governing it. We will then find 
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analytic expressions for the steady-state current and the 
rate of heat produced as a function of the model param- 
eters. In section III, we will consider our model in the 
presence of an external load and find the steady-state 
current and the rate of heat flows. In section IV, ex- 
pressions for the efficiency and coefficient of performance 
(COP) will be determined depending on whether the en- 
gine works as a heat engine or as a refrigerator. Regions 
in the parameter space where the model works as a heat 
engine, as a refrigerator and as neither of the two will be 
determined. We will explore how current, efficiency and 
COP behave as the model parameters vary. We will also 
compare the efficiency of the engine at two operating con- 
ditions (at maximum power and at optimum value with 
respect to energy and time) with those of endoreversible 
and Carnot engines. Lastly, we summarize and conclude 
in section V. 



II. ZERO EXTERNAL LOAD 

The model we take is a modified version of the one con- 
sidered by Jarzynski and Mazonka 5] in modelling Fcyn- 
mann's ratchet and pawl system. Our model is meant 
to capture the motion of a particle in a ratchet potential 
due to the thermal kick it gets from periodically placed 
hot and cold reservoirs along the path. 

Consider a particle that moves by hopping on a one- 
dimensional lattice, with lattice spacing d, assisted by 
the thermal kick it gets from periodically placed hot and 
cold reservoirs along its path. The particle is also exposed 
to an external discrete sawtooth potential which has the 
same period as that of the reservoirs. In one cycle, the 
particle walks a net displacement of three lattice sites, 
either to the right or to the left. This corresponds to the 
particle crossing one sawtooth of a sawtooth potential. 
The potential energy at site i, Ui, where i is an integer 
that runs from — oo to +oo, is given by 

U l = £[z(mod)3- 1]. (1) 
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FIG. 1: Plot of discrete sawtooth potential without load. 
Sites with dark circles are coupled to the hot reservoir (Th) 
while sites with open circles are coupled to the cold reservoir 
(T c ). Site 1 is labeled explicitly and d is the lattice spacing. 



Here E is positive constant having a unit of energy. 
The temperature profile at site i, Tj, is given by 
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where Th and T c are the temperatures of the hot and 
cold reservoirs, respectively. Figure 1 shows the values of 
sawtooth potential and of the temperature at the given 
lattice sites. 

The jump of the particle from one lattice site to next 
lattice site is assumed to be random in nature. The 
jump probability is determined by the amount of energy 
it crosses and the temperature of the heat reservoir to 
which it is coupled. Accordingly, the jump probability 
per unit time of the particle making a jump from site i 

-AE 

to site i + 1 is given by Te T * , where AE = Ui+i — IJi 
and r is the probability that the particle will attempt a 
jump per unit time. We take Boltzmann's constant, kg, 
to be unity. When the particle attempts to jump, first 
it decides which way to jump (either to the left or right) 
with equal probability and then jumps according to the 
Metropolis algorithm if the value of AE < 0, then 
the jump definitely takes place; if AE > then the jump 
takes place with probability exp( ~^ E ). 

The dynamics of the model can be studied by mapping 
the model to a spin-1 particle system which exhibits iden- 
tical behavior 0. Denoting the spin states, s, by (0, ±1), 
the energy function of the spin-1 particle is defined by 



E(s) = Es. 



(3) 



The change of state of spin s corresponds to the jump 
of the particle. If s changes from -1 to or from 
to 1 or from 1 to -1, then this is the same as the 
particle jumping to the right. The reverse process corre- 
sponds to jump to the left. The three possible states and 
their corresponding energy values is shown in the Table I. 

The dynamics of the particle is then described by 
stochastic jumps among the three states. The process 
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TABLE I: The three states of the system and their corre- 
sponding energies 



is Markovian and we can describe the evolution of the 
states with rate equations. Let the probability for the 
system to be found in state n at time t be given by p n (t). 
The rate equations governing the evolution of the three 
states are 



~Jjp = (Pnn'Pn' ~ P n 'nPn),n, n' = 1, 2, 3. (4) 
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P n ' n is the transition probability rate at which the sys- 
tem, originally in state n, makes transition to state n'. 
Here P n > n is given by the Metropolis rule. For example, 
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(5) 



In the above expressions, the factor | is due to the de- 
cision for the particle to jump either to the left or to 
the right. The rate equation for the model can then be 
expressed as a matrix equation 
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where p = (pi,P2,P3) T ■ Here, R is a 3 by 3 matrix which 
is given by 
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where fi = e T <= and v — e T <» . Note that the sum of 
each column of the matrix R is zero, R ran = 0. This 
shows that the total probability is conserved: Yl n P n ~ 

£(i T .p) = i T .(rRp) = o. 

The steady state probability distribution p of Eq. (4) 
is obtained by solving Hp = 0. We find the normalized p 
to have components given by 
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The presence of the hot and cold regions along the 
lattice leads to unidirectional steady state current, J. 
This steady state current, J, can be found as the dif- 
ference between the current towards the right, J + , and 
the current towards the left, J~ , between any two states: 
J = J + — J~ . Selecting processes taking place between 
states 2 and 3, the current towards the right, J + , is given 
by 



T(i?32P2), 



(11) 



while the current towards the left, J , is given by 
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(12) 



After some algebra the explicit expression for the current 
J takes a simple form 



2(2 + +M + /J 2 )' 



(13) 



In each cycle, the particle walks a net displacement of 
three lattice sites, 3d. Therefore the drift velocity, v, of 
the particle is 



v — 3dJ. 



(14) 



Notice that the net current is to the right as long as 
Th > T c and zero when Th = T c . 

Let us next find the amount of heat transfer per cycle be- 
tween the hot and cold reservoirs as the particle climbs 
up or down the potential. We assume the case where 
there is no energy transfer via kinetic energy due to par- 
ticle recrossing of the boundary between the hot and cold 
reservoirs 7, 8J. When the particle jumps from state 2 to 
state 3, it takes heat from the hot reservoir whose amount 
is sufficient to climb up the potential energy difference 
between the states and equal to E. When the particle 
jumps from state 3 to state 2, it gives heat to the hot 
reservoir by losing its potential energy and is equal to 
E. Thus, the net heat per unit time taken from the hot 
reservoir due to climbing up or down the potential, Qh, 
is given by 



Qh = ET(R 32 p 2 ~ R23P3)- 



(15) 



After substituting the values of p 2 and p 3 from Eqs. 
(9)and (10) in Eq. (15), we get 



= ET 



(2(1 + M + M 2 ) (2 + *,))' 



(16) 



When the particle jumps from state 3 to state 1 and 
from state 2 to state 1, it gives heat to the cold reservoir. 
When it jumps from state 1 to state 3 and from state 
1 to state 2, it takes heat from the cold reservoir. The 
net heat per unit time given to the cold reservoir due to 
climbing up or down the potential is given by 

Q c = ET(2R 13 p 3 - 2R 31 p 1 + R 12 p 2 - R 21 p 1 ). (17) 



FIG. 2: Plot of discrete sawtooth potential with load along 
with the discrete temperature profile on the lattice 



After substituting the values of p^'s, we obtain 

mO - m) 



= ET 
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The above results show that Q c = Qh- This clearly shows 
that a heat taken from a hot reservoir directly goes to 
the cold reservoir without doing any work. Unlike this 
work, in our previous work £| , there is work done against 
viscous medium even in the absence of external load. 
Let us now compute the rate of entropy production. The 
rate of entropy production related with the flow of heat 

from the hot reservoir is given by Sh = -jr^ while the 
rate of entropy production related with the flow of heat 

to the cold reservoir is: S c = The net rate of entropy 

production given by S — Sh + S c takes of the form, 
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Notice that S is non-negative which implies that the sys- 
tem is consistent with the second law of thermodynamics. 



III. 



NON-ZERO EXTERNAL LOAD 



Let us consider the model in the presence of constant 
external load, /, added to the sawtooth potential as 
shown in Fig. 2. Accordingly, the potential energy will 
now be changed from Ui to Ui + ifd. The method of 
solving for steady state behavior is the same as that for 
the load-free case. We consider our model when the load, 
/, is greater than zero and for Metropolis algorithm to 
hold, we need to limit the range of / to < / < . For 
this range of the load, we find the matrix R, 
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where a — e T " and b = e T '» . Note that the sum 
of each column of the matrix R is zero, which shows 
that the total probability is conserved. The steady state 
probability p satisfies the matrix equation, Hp = 0. We 
solve for p and after normalization the final results are of 
the form: 
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Using similar approach as in section II, the expressions 
for the steady state current, J, and the drift velocity, v, 
are respectively found to be given by 
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When the Brownian particle walks along the potential 
with additional load, it takes heat from the hot reservoir 
and gives some part of it to the cold reservoir and uses 
the rest for climbing up the load. The difference between 
the rate of heat energy that the Brownian particle takes 
from the hot reservoir, Qh, and the rate of heat energy 
that it gives to the cold reservoir, Q c , is the rate of useful 
work, W, that the particle uses to lift the load; i.e, 



W = Q h - Qc = fv. 
From Eqs. (25) and (26), we get 



W= fv = 3fd 
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The rate of heat energy taken from the hot reservoir by 
the particle, Qh, can be obtained by a similar approach 
as in section II and we find its expression to be 



Q h = (E + fd) 
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The rate of heat energy given to the cold reservoir ■ 
given by 



Q c = (E-2fd)- 



Tn(bav - f ) 



2(2 + i/6)(l + a^+^) 

The difference between Qh — Qc will be 

T^{bav - H) 

Qh - Qc = 3fd— — ~ T^T) 
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which is exactly equal to W as given in Eq. (27). 
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FIG. 3: Plot of j versus A for r = 1 and e — 2 



IV. THE MODEL AS A HEAT ENGINE, AS A 
REFRIGERATOR AND AS NEITHER OF THE 
TWO 

To specify the model, one requires to specify the six 
quantities: T, d, E, T c , Th and /. Taking T, d and T c 
fixed we still have three parameters E, Th and / that can 
be varied independently. We convert these into three 
dimenssionless parameters e, r and A where e = M-, 

+ c 



In addition, we introduce a 



t = | - 1 and A = |. 

dimenssionless current j = p. We take > T c for the 
rest of our work. 

The current j is now a function of e, r and A. Figure 
3 shows a plot of j versus A for fixed values of e and r. 
The figure shows that the current is positive as long as the 
load is less than a certain value of A . This corresponds 
to the region where the model works as a heat engine. 
Using Eqs. (28), (29) and (30) the efficiency of the heat 
engine, 77, takes the expression 
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(31) 



On the other hand, when the load is large enough the cur- 
rent becomes negative and that implies that the model 
works as a refrigerator. The COP, P re f, of the refrigera- 
tor then takes the expression 
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(32) 



From this equation, Eq. (32), we note that in order for 
the model to function as a refrigerator the upper limit 
for A must be |. 

The set of points in the parameter space at which cur- 
rent changes its direction differentiates the domain of op- 
eration of the model as a refrigerator from that as a heat 
engine. Using the expression for J, Eq. (24), the value 
of A at which the current reversal takes place is given by 



A = 



(2t + 3) ' 



(33) 



Note that this value of the load where current is zero 
is usually called the stall force for molecular engines 



FIG. 4: Plots showing the three regions of operation of the 
model in the (a) A-r paramters space for e — 4 and in the (b) 
A-e paramters space for r = 1 

0- When we evaluate the values of both 77 and P re f 
as we approach this boundary determined by Eq. (33), 
we find that they are exactly equal to the respective 
values for the Carnot efficiency and the Carnot COP: 

' and lim P re f ~~ 



lim 77 



...... lim ^ , , , 

T h J^o- ' LJ (T h -T c , 

clearly demonstrates that the boundary at which current 
is zero corresponds to the quasistatic limit be it from the 
heat engine side or from the refrigerator side. 

Let us introduce dimenssionless parameters qh = 

and q c = -py-. In analyizing the operation of our model 
above, we have identified that when 
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the model works as a heat engine while it works as a 
refrigerator when 
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Therefore, the model neither works as a heat engine nor 
as a refrigerator when 



A> 2- 



(36) 



Figure 4a shows the three regions in the A-r parameters 
space in which the model operates as a heat engine, a 
refrigerator and as neither of the two for fixed value of 
e = 4. On the other hand, Fig. 4b shows these three 
regions in the A-e parameters space for fixed value of 

T= 1. 

Let us now further investigate how the current, efficiency 
and performance of the refrigerator behave as a function 
of the different parameters characterizing the model. The 
plot of 77 versus A shows that, the efficiency , 77, increases 
with increase in A until it attains its maximum value 
(Carnot efficiency) as shown in Fig. 5a. On other hand, 
the plot of P re f versus A shows that within the range 
where the model works as a refrigerator, coefficient of 
performance of the refrigerator, P re f, decreases from its 
maximum value (Carnot refrigerator) as A increases (see 
Fig. 5b). 

Figure 6 shows how the current j behaves as a function of 
e. The figure shows the presence of maximum current at 



FIG. 5: (a) Plot of 77 versus A for e = 2 and r = 1 (b) Plot 
of P re f versus A for e = 2 and r = 1 
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FIG. 6: Plot of j versus e for A = .2 and r = 1 



a certain finite value of e for a fixed r and A. This point 
corresponds to maximum power delivery at which the en- 
gine operates with maximum power efficiency, t\mp- Let 
us compare the value ijmp of our engine with the corre- 
sponding value of an endoreversible engine. For indorc- 
versible engine which exchanges heat linearly at a finite 
rate with two reservoirs, Curzon and Ahlborn [Tol Hi) 
showed that the efficiency at maximum power, t]ca> is 
given by 



VGA 




(37) 



Figure 7 shows how the efficiency of our model heat en- 
gine when it operates with maximum power, t]mp, and 
that of t\ca behave as a function of r. The plots show 
that as t increases the gap between t]mp and r\cA de- 
crease and coincide at a finite value of r. Then their 
difference get larger as r increases. When we compare 
the two efficiencies, t\mp and 77c ^4, f]CA is found by as- 
suming linear heat conductivity while t]mp is obtained 
numerically without taking any assumption. This illus- 
trates that rjMP works for the entire range in the allowed 
parameter space while r\cA is specific and has limited 
significance. 

Let us next compare Carnot efficiency, tjcar, with that 
of optimized efficiency, tjopt- The 770 pt for our model 
can be found using the argument stated by Hernandez. 
et al 03. We briefly summarized the method Hernandez. 
et al |l2j | in our earlier work |4j . The optimized efficiency, 
Vopt, lies between maximum efficiency and efficiency un- 
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FIG. 7: Plots of 7]ca, tjmp, tjopt and t\car versus r, where 
the model engine is put to function at A = 0.2 while e is fixed 
depending on whether it is working at either maximum power 
or optimized efficiency 
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FIG. 8: Plot of Q versus e for A = .2 and r = 1 

der maximum power and it is given by optimizing 

SI = 2W- ^ Th ~ Tc) Q h . (38) 
J-h 

The plot of Cl versus e in Fig. 8 shows that the function 
definitely has optimum value at finite e. Evaluating the 
efficiency at this particular point in the parameter space 
gives us the optimized efficiency, tjopt- We plot tjopt 
and rjcAR in the same Figure 7 that we plot the other 
efficiencies. The plots of t\opt and t\car versus r show 



that rjoPT lies between tjmp and t\car- This undeniably 
illustrates that the operation of the engine at optimized 
efficiency is a compromise between fast transport and 
energy cost. 



V. SUMMARY AND CONCLUSION 

In this work, we introduced an exactly solvable model 
of a heat engine with minimum ingredients. We obtained 
closed form expressions for current, efficiency and COP 
of the model. We showed that at quasistatic limit the 
values of both efficiency and COP go to that of Carnot 
efficiency and Carnot COP, respectively. We then ex- 
plored the basic properties of the microscopic heat engine 
by varying the parameters describing the model. We fur- 
ther studied the efficiency when the engine operates with 
maximum power and compared this efficiency with those 
values one gets by using the finite-rate linear heat ex- 
change assumption of Curzon and Ahlborn jToL ITlj . The 
results of optimized efficiencies of the model are also re- 
ported. It is worth to note that the particle walks in 
non-viscous medium. Hence the model does not work 
as a heat engine in the absence of external load. In the 
presence of external load, even though the model and its 
corresponding dynamics is completely different from the 
one we studied earlier, there is a qualitative agreement 
between this work and the earlier work. Therefore, one 
can take this work as independent check of the results we 
found in the previous work on Brownian heat engine Q . 
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